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1. Introduction and Summary 

Since time immemorial [1] people have looked for liquid crystal phases with bond- 
orientational order [2] . Aside from the possibility of a new type of liquid crystal, hexatic- 
type order can be an intermediate stage in the continuous freezing of a liquid in three 
dimensions. Much more recently there has been an explosive progress in new chiral phases 
in both two [3,4,5], and three dimensions [6,7]. In this paper I propose both a uniform and 
defect laden ground state of a liquid crystal with both hexatic order and chirality. 

Toner [1] has proposed that nematic liquid crystals, upon cooling, could form a liquid 
crystalline phase with nematic order and hexatic order in the plane perpendicular to it. 
In section 2 I will consider the melting of a chiral columnar phase [7] into a chiral liquid 
with hexatic order. I discuss both the Landau theory as well as the melting of the crystal 
through dislocation loop unbinding. When a proliferation of edge and screw dislocations 
develops, the crystal melts, leaving a normal liquid with hexatic order and no residual 
constraints associated with the crystal rigidity. Typically, long thin molecules do not have 
columnar phases, but freeze via smectic phases. I expect then that chiral, disc shaped 
molecules could participate in the structures discussed here. 

It is well known [8] that cholesteric phases, under applied magnetic fields can unwind 
as well as tip up and become conical phases. However, the conical phases exist only when 
K 2 > K 3 , which is seldom the case experimentally. When K 2 < K 3 there is a discontinuous 
transition from a cholesteric with director in the xy-plane to a nematic aligned along the 
field in the z direction. In section 3 I find that hexatic order has the same effect as the 
magnetic field, although it does not participate in unwinding the helical pitch. Adapting 
the known consequences of fluctuations and nonlinearities of polymer nematics [9] and 
cholesterics [10] to the chiral polymer hexatic, I find anomalously large elastic constants 
which depend on the polymer length. For sufficiently long polymers K 2 will eventually 
grow larger than K%. If liquid crystal phases can still exist for these long polymers, Meyer's 
conical phase would appear. This would be evidence of hexatic order. When the nematic 
director has a conical texture, as in a smectic-C* texture, and the molecules are chiral, 
an electric polarization can develop. Thus this phase could exhibit ferroelectric liquid 
crystalline behavior [11] . 

Since the bond order is defined in the plane perpendicular to the nematic director 
n, its value must be determined in a way which takes into account the texture of the 
director. Because there is a nematic director n and a bond-angle order parameter 9q, the 
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non-chiral theory is similar to that of the A phase of 3 He. In section 4 I will discuss the 
geometric constraints relating the nematic director to the bond-order field in the absence 
of free disclinations. This is known as the Mermin-Ho relation in 3 He [12]. It plays an 
important role in understanding the allowed ground states of the chiral N+6 phase and 
elucidates the structure of the allowed defects. 

With chirality added the N+6 order can exhibit a large variety of phases similar in 
structure to known phases of smectics. When chiral molecules participate in N+6 order, 
two kinds of twisting are allowed [7]: twisting of the nematic order, leading to cholesteric 
states, and twisting of the bond order, leading to braided moire states. The twisting bond 
order, with an ever increasing angle, is analogous to the layered order of smectics with 
the pitch of the bond order equivalent to the layer spacing of the smectic. This allows me 
to consider defect phases similar to the Renn-Lubensky TGB state in chiral smectics. I 
discuss this in section 5. 

2. Landau Theory and Melting of the Hexagonal Columnar Phase 

I first derive the free energy of the chiral N+6 phase as the sum of its constituent 
parts. Fluctuations in the nematic director are described by the Frank free energy density: 

f fl = |(V-n) 2 + ^[n.(Vxn)- go f + |[nx(Vxn)f (2.1) 

where Ki are the Frank elastic constants and 27r/<zo is the equilibrium cholesteric pitch. 
In order to define bond-angle order, I must, everywhere in space, define a right-handed, 
orthonormal triad {ei,e2,n}, where e.2 = n x e±. In this case, the usual hexatic order 
parameter, ifj^ = \ipe\e 6 is defined by 

iSP(n,r) 

where 9i(r) is the angle between the particle % and the basis vector ei as measured around 
r with the sign of the angle determined by n, and -P(n, r) is the set of particles that are in 
the plane perpendicular to n through the point r. The vector ei is chosen unambiguously 
throughout space (possibly modulo rotations by 27r/6) and may be thought of as pointing 
to one of the ground state nearest neighbors. More precisely, 



and the sum (2.2) is taken only over particles at position i-j such that n(r) • (r; — r) = 0. 
The bond-angle order parameter now depends on the texture of n. As a result, one must 
take into account the nematic when taking derivatives of 9$: a covariant derivative must 
be employed. For now I will ignore this complication and will return to it in Section 4 and 
will justify the naive analysis a posteriori. 

Because the definition of 9q depends on the direction of n, under the transformation 
n — > — n, (2.3) implies that Qi — > — 0, and hence 9$ — > —9q (or equivalent ly, tpQ <-> ijjg). 
When constructing a free energy that includes 6>6, the overall nematic symmetry n — > — n 
must be preserved. Owing to the definition of 9q, the free energy can include any term with 
even powers of n and 9§ together. The spin wave theory for 9q includes spin-stiffnesses as 
well as a new chiral term: 

= K *~ 2 K * (n • W 6 ) 2 + (W 6 ) 2 - K%n • W 6 (2.4) 

where, since there is a preferred direction n, I have included the possibility of anisotropic 
stiffnesses. The final term is chiral and indicates the tendency for the bond-angle order 
parameter to rotate around the nematic director with pitch 2-7r/go [7,13]. Finally there are 
additional non-chiral couplings between n and 9q [1,14]: 

F A e 6 =C(h- V6> 6 ) [n • (V x n)] + C'V9 6 • V x n (2.5) 

The free energy is the sum of terms F = J d s x {jF ft + J r d6 + jF fi g 6 }. 

Recently the theory of chiral polymer crystals has been formulated [7]. One would 
imagine that the total free energy resulted from the melting of such a crystal. Indeed, this 
is the case. Considering only quadratic fluctuations around a ground state with (ft) = z, 
the additional free energy coming from the two-dimensional crystal displacement field u is 

^"crystal = K d i U j ~ e *A) 2 + + I*' (d x Ui ~ SUi) 2 (2.6) 

At long wavelengths 9q and 5rii are locked into crystal deformations by 9q = ^eijdiUj and 
Srii = d z Ui. Thus 28 z 9q — eijdiSrij = 0. When dislocations are introduced into the crystal, 
derivatives do not commute so d z di ^ did z . In this case I introduce w^i which is equal to 
dyUi away from dislocations. As in [15] I introduce a dislocation density tensor a^, the 
density of dislocations running in the fi direction with Burgers vector in the % direction 
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and e^^d^Wji = —a^. With drii = w z i and 9q = \tijWij the long wavelength constraint 
between 6 and 5rii becomes instead [7] 



2<9 2 6> 6 - tijdibnj = -Tr[ct] 



(2.7) 



where eij is the two-dimensional anti-symmetric tensor. In particular (2.7) constrains the 
two chiral terms to differ only by the density of screw dislocations. 

Since a proliferation of edge dislocations running along the z direction would be suf- 
ficient to melt a two-dimensional columnar crystal, one might think that the constraint 
would not be altered and that the free energy of the liquid state would be augmented by it. 
However, this is not the case. In the usual scenario of two-dimensional melting [2] dislo- 
cation pairs unbind, leading to a hexatic liquid. In the case of a columnar crystal, parallel 
dislocation lines can unbind, leading to a hexatic. In a columnar crystal, defects running 
parallel to the columnar direction (z) are edge dislocations. Those defects lying in the 
perpendicular plane (xy) are either edge or screw dislocations, depending on whether the 
Burgers vector is normal or parallel to the dislocation line, respectively. Thus dislocation 
loops which contribute to the melting of the crystal through edge dislocations along z, will 
contain, as well, screw-like and edge-like dislocations in the plane. Since the energy per 
unit length of a pair of edge dislocations parallel to z is finite, infinite dislocations, ending 
at the boundaries, will not proliferate. Instead, dislocation loops will unbind [15] leading 
to a proliferation of screw dislocations. To incorporate the presence of screw dislocations 
I add to (2.6) a term representing the dislocation free energy density [16]: 



where represents the core and elastic energy of the defects. In the pure crystal 

the energy per unit length of an unpaired edge dislocation parallel to the i-axis diverges 
logarithmically with the system size. Dislocation loops, however will have a finite energy, 
scaling like RlnR where R is the size of the loop. Below the melting transition, where the 
defect loops are bound, the free energy will favor a^i = and hence Tr[a] =0 everywhere. 
Thus 2<9 Z 6> 6 = V x 5n. However, when dislocation loops unbind, the effective quadratic 
energy for ot^i will be finite. The diagonal part of (2.8) becomes 



Thus the constraint (2.7) will be demoted to merely a preference for configurations with 
2<9 2 6>6 = eijdidrij [17,15] and can be absorbed into shifts in k\, K 2 and C. Thus the 
melted crystal is precisely described by our liquid crystal free energy. 



bT = S M i^(q)a M i(q)a^(q) 



(2.8) 



fc^diag = ^screw {Tr["]} 2 = ^screw [ld z & - V X 5n] 2 



(2.9) 
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3. The Uniform Conical Phase 

I now look for ground states as a function of Ka, the amplitude of the hexatic order 
parameter. The results of this section closely follow those found by Meyer for a cholesteric 
in a magnetic field parallel to the pitch axis [8]. The analysis here, however, involves the 
bond order parameter 9q and its equilibrium configuration as well as that of the nematic 
director. For simplicity I take = K A = Ka and likewise take C = 0. I consider as a 
class of ground states the conical states proposed by Meyer [8] : 

h = [cos 4> cos qz, cos 4> sin qz, sin 0] . (3.1) 

where q and (p are free constant parameters. The equations of motion for 9q are 

-K A V 2 6 6 - K A q V • n = (3.2) 

and so V 2 6> 6 = as n is divergence- free. Thus the only solutions for 6> 6 are linear functions 
of the coordinates and so v = V6>6 is a constant vector. Inserting the ansatz (3.1) and a 
constant vector v into (2.1), (2.4) and (2.5), I have 

F = n{ ^v 2 - K A q sin 0v o • z + ^ [q cos 2 - q ] 2 + sin 2 cos 2 j (3.3) 

where O is the volume of space and the oscillating terms drop out upon integration over 
space. Minimizing with respect to v I find 

v = V0 6 = ^ o sin0 (3.4) 

and thus, with x = cos 2 4> 

F eS = n |^ (qx - q f + ^q 2 x{l - x) - ^q 2 (l - *)} (3.5) 
This is precisely the energy studied by Meyer for a cholesteric in a magnetic field of 



strength H = (Jq\/Ka- In [8] the possibility of the cholesteric unwinding by changing 
q was considered along with the conical state. In this case the former effect can never 
happen: because v must be constant, the effective field will always be parallel to the 
pitch axis. Translating Meyer's results I find that for K 2 < K3 there is a transition with 
increasing v^AQo from the pure cholesteric state ((f) = 0, q = qo) to a pure nematic state 
((J) = 7r/2) at \fK~Aq G = q^^fK^. For K 2 > K 3 there is an intermediate phase as well. 
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For y/TZXqo < q y/K^ the cholesteric state persists, while for y/K^qo > %K 2 /\/K^ the 
nematic state persists. In between 

q = qoVK A /K 3 (3.6) 

and 

cos2 ^ i= WTO(*ZM^ (3 .7) 



for qo\/Ks < \fK~Aq~o < qoK 2 / 'y/K^. Thus the conical phase should appear continuously 
as a function of q~o, or, alternatively, as a function of for fixed %. Hence as the hexatic 
order grows the cholesteric state will become conical and, eventually, nematic. The conical 
state has, in addition to a smectic-C*-like director texture a rotating bond order as well, 
with pitch q . If qo/Qo is n °t rational, the fluctuations around this ground state should be 
akin to those in an incommensurate smectic [16,18]. 

It is unlikely that q~o = as it is allowed by the same symmetry that allows the 
cholesteric coupling q . However, even if this is the case I expect C ^ 0. In the case 
q = 0, the C acts as a generated q , in some sense. I will still have a magnetic field type 
term, though now it depends on q and 0, which should produce a conical phase. The phase 
diagram at qo = will depend on all of K 2 , K3, Ka and C, and mapping it out would be 
straightforward yet tedious. 

Turning back to q 7^ 0, we have found that the conical phase can persist for K 2 > 
K 3 which is usually not the experimental case. However, in light of recent work on the 
fluctuation enhancement of elastic constants [9,10] one might hope that this inequality 
could always be met. Indeed since there is no constraint locking 2d z 9§ with x 5n the 
non-chiral part of the free energy from [10] would be 



1 I 2 
V x -it- -(d z u) 2 +u 



+ Gp 2 (d z u) 2 +^[d z u] } (3.8) 



where E is the bulk compression modulus, po the average polymer density, d z u = Sn, 
G = kT£/po, I is the typical polymer length and J 7 a is expanded to quadratic order in 
8n where n rs z + 5n. It was found in [10] that for polymers longer than £q, E, K 2 
and K3 took on anomalous, ^-dependent values. For polymers with only steric entropic 
interactions £q = Lip /a 2 where Lp is the polymer persistence length and a is the average 
areal interpolymer spacing. It was found that K 2 (£) ~ (£/£o) - 20 and K 3 (£) ~ (£/£ ) 0A5 
to two loops in a d = 4 — e expansion at d = 3. 
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To the free energy in (3.8) I add the hexatic stiffness J r o 6 . Expanding around n = z I 
can consider the possible corrections to the hexatic stiffness Ka as well as the chiral cou- 
plings qo and go ■ Following the analysis in [10], I find that Ka and go do not anomalously 
renormalize and suffer only finite shifts in their values. To illustrate, I consider adding 
Tq q to (3.8) with K\~ = k\ = Ka- I am interested in fluctuation-induced changes to 
Ka and qo. As in the argument of [10] for the non-renormalization of go, I consider the 
go ~ limit. Thus I do not consider corrections to quadratic propagators arising from the 
coupling KAqod z iI ■ V ± 9q. It is easy to see that, since I consider only corrections linear 
in go, Ka does not acquire any infinite renormalization, and, in fact, even if I consider 
K^ 7^ k\ the corrections from fluctuations in 5n only lead to finite renormalizations. 
The corrections to go must come from non-linearities in the chiral term itself; no other 
term can generate the appropriate chiral coupling. The tadpole graph arising from this 
is proportional to (5n 2 (r)) and thus only contributes to a finite shift in Ka% in three 
dimensions, where nematic order persists. Alternatively, by a judicious choice of rescal- 
ings [19], one can show that the longitudinal part of u is irrelevant to the renormalization 
of the Frank constants, and thus the new coupling, which only involves the longitudinal 
part of -u, will not change the results in [10]. Thus the mean field analysis of the uniform 
conical state holds when I replace the constant Frank constants with ^-dependent Frank 
constants. Hence, for sufficiently long polymers, it may be possible to have the unusual 
situation of K 2 > K 3 , though since K 2 /K 3 ~ (£/£ ) - 05 , a twofold increase in the ratio 
of Frank constants implies the polymer must be 10 6 £q. For instance, in DNA Lp = 600A 
and a = 35A lead to £o w 18/um. With typical chain lengths on the order of centimeters, 
this can lead to a 37% increase in the ratio of K 2 to K 3 . This increase leads one to hope 
that the conical phase could persist and could be observed. 

Aside from the light scattering behavior of this phase due to its N+6 character [1], it 
should be possible to see this phase through crossed polarizers. As hexatic order grows past 
the cholesteric-to-conical point, the extinction directions (perpendicular to the director) 
will continuously change from the nematic direction to the (perpendicular) cholesteric 
plane. This will happen without a density modulation forming, as in the smectic-C* phase 
which would have a similar behavior. The density modulation could be detected (or not) 
via electron microscopy and x-ray scattering, thus distinguishing this new phase from the 
smectic. As in the moire phase of chiral columnar crystals [7], in each constant z cross 
section the structure function will be, in Fourier space, that of a hexatic - six broad spots. 
When looking at the structure of the bulk sample, however, the rotation of the bond order 
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will merge these spots into a ring in the q x -q y plane. The periodicity of the rotation will 
lead to additional rings at q z = 6nqo sin where n = ±1, ±2, .... 

In addition, as in a smectic-C*, this uniform conical texture can exhibit ferroelectric 
behavior through a term in the free energy 

JF^ = n (v • n) E • (v x n) = /jq sin cos <j> [E y cos — E x sin gz] (3.9) 

This extra interaction preserves the nematic symmetry and is chiral, since under parity, 
E — > — E. I have taken the "layer normal" to be v in constructing (3.9). There is, however, 
a term that is allowed in these chiral phases that has no analogue in the smectic-C* . Since 
the direction of v is determined by the chiral coupling q and n, a non-chiral term is 
allowed 

J^flc = pE ■ (v x n) (3.10) 

which preserves the nematic symmetry (n and v both change sign) and is not chiral. Under 
spatial reflection v x n unambiguously changes sign, as does E. This term is allowed 
because there is a difference, for fixed n, between a configuration with v and —v. They 
differ by the handedness of the hexatic helix. In a smectic-C* the layer normal does not 
have an unambiguously defined direction because of the up-down symmetry of the layered 
structure. Unlike a smectic, this ferroelectric liquid crystal does not have a density wave 
as in the usual ferroelectric smectic liquid crystals. The wave that supports the conical 
texture is one of twisting hexatic order along v. In section 5 I will make the analogy 
between smectics and chiral hexatics more precise. Note that jF flc also implies that an 
applied electric field will favor configurations in which the hexatic order twists in planes 
perpendicular to n. 

4. Topological Constraints: The Mermin-Ho Relation and Saddle-Splay 

Up to this point I have ignored the topological relation between 6 and n. From the 
definition of 6> 6 arising from (2.3), it is clear that textures in n can influence the local value 
of 9q. The two orientational order parameters are related to each other in much the same 
way that the two order parameters of 3 He-A are related, namely, through the Mermin-Ho 
relation [12]. I will sketch here a brief derivation, following [20]. Let {e!,e 2 ,n} be a 
right-handed, orthonormal triad. If I consider a rotation 59q about n, then 

e[ + ie' 2 = e lS0e (e x + ie 2 ) (4.1) 
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and so 



e'i -ei « -59 6 e 2 
e 2 — e 2 ~ ^^ei 



(4.2) 



Thus 56> 6 « e! • 5e 2 , or 



(4.3) 



Taking the curl of both sides and using e! x e 2 = n, we have 



[V x v] 



1 



(4.4) 



the celebrated Mermin-Ho relation. The content of this relation is that in the ground state 
of the N+6 phase, disclinations in 6> 6 are forced in by the nematic texture. It quantifies a 
simple topological consequence of non-cylindrically-symmetric molecules. For concreteness, 
take a biaxial molecule in a cylindrical geometry with normal boundary conditions. One 
might imagine a +1 disclination line will run up the center of the cylinder. However, 
it is well known that such a geometry will "escape into the third dimension" along the 
disclination. However, if the molecules have the same (biaxial) orientation at the cylinder 
wall, when they escape they will be compelled to force a disclination into the biaxial order 
parameter in the center. This is illustrated in Figure 1. This frustration is summed up in 
(4.4). Note that in the previous section, I took V x v = 0, which, in the case of the Meyer 
ansatz, is justified by simply plugging (3.1) into (4.4). 

When I look for ground state configurations I must take the constraint into account. 
The excited states differ from the ground state by the energy of a free disclination. In 
the hexatic phase, this energy diverges with the system size, so the Mermin-Ho relation 
should be enforced. This is not unlike hexatic phases on flexible membranes. There 
disclinations were forced in by gaussian curvature [21]. In particular, if I take as the the 
hexatic director N = cos 06©i + sin 6>6e 2 and compute the covariant derivatives of N in the 
local basis determined by {e l5 e 2 , n}, I find 



^ = — [D t N 3 f + K - [D, Nj f = K — (W 6 - O) 2 + ?—JL- [n ■ (W 6 - Q)] 2 (4.6) 



D^Ni = -e„ ■ [V6> 6 -Q] sin ^ 6 
D^N 2 = V [W 6 - ft] cos # 6 



(4.5) 



where n = £3 and O 



= e^d^e^- The appropriate elastic free energy would be 
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In terms of N, the term n • V6> 6 , for instance, becomes n^jNiD^Nj = n • (V6> 6 — Q). 
Thus, in general, one should always include with derivatives of 9§ as a consequence of 
the wobbling about of n. Because a change of basis vectors ei(r) and e 2 (r) should not 
change the physics of 6>6, I will consider only "gauge-invariant" correlations of 9q through 

G M6 (r; r) = ( e 6 ^« exp{z J fi ■ M}e~ 6ie ^ ) (4.7) 

where T is a path connecting r and 0. When this Greens function is non- vanishing for large 
separations r, there is long-range hexatic order. This correlation function will be invariant 
under shifts of O by Vuj and 6> 6 by uj, where uj is an arbitrary function of space. These 
transformations are precisely those which do not change VxO and hence do not change 
the disclinations structure of 6> 6 . When V x = everywhere, there are no obstructions 
to choosing e;(r) to be constant in space, leading to O = and to the usual definition 
and interpretation of 9$. Note that when the nematic is highly disordered, with a rich 
wobbling texture, the contribution to (4.7) from the nematic curvature can greatly reduce, 
and possibly destroy, long-range hexatic order. 

Since unbound disclinations cost a logarithmically-divergent energy they will not ap- 
pear in equilibrium. However regions where Vxfi^O have precisely the same divergent 
energy as a free disclination. Writing 9$ = #j; mooth + 9g ing we can expand around a back- 
ground of disclinations chosen to screen these regions with V#g mg = O, leaving only smooth 
variations in 9q. The dislocations are locked to Vxf! according to the Mermin-Ho relation 
(4.4). Notice that in the case that n corresponds to a field of normal vectors to surfaces, 
n(r) • [V x O(r)] is the Gaussian curvature of the surface at r to which n(r) is normal. 
Rewriting the curvature, I find that n ■ (V x fl) = —V ■ [(n • V) n — n (V • n)], the well 
known saddle-splay term of nematic liquid crystals. Thus I can identify C' with X24 the 
saddle-splay elastic constant. 

This can be understood geometrically: to lowest order in 8n 

[V x V9 6 ]i = ^eijke^n^djn^dkn 1 m eij k djbn x d k bn y (4.8) 

Disclinations pointing along z are constrained by 

[V x V6> 6 ] z = d x Sn x d y Sn y - d x Sn y d y Sn x (4.9) 

in other words, disclinations must appear if the directors are the normals to saddle-surface. 
The saddle-splay term measures the amount of saddle-like deformations in the director. 
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Since 6> 6 is single-valued on a surface, the Mermin-Ho relation tells us that the saddles 
must be such that the integral of their saddle-splay around their boundary is a multiple 
of 27r/6. However, in the N+6 case, the director n need not be the collection of normals 
to a surface. Indeed, in the double twist texture of liquid crystal blue phases there is no 
surface to which the directors are normal. If such a texture were present in the chiral N+6 
system, (4.9) implies that as a circuit is traversed around a central molecule of a double 
twist cylinder, there would be fewer than z = 6 nearest neighbors since the saddle-splay is 
positive. Considering the director texture around the center of a double twist core, Figure 
2 shows 19 braided polymers with the center one being straight. It is not necessary that 
these be polymers, just that the nematic director be tangent to the lines. This texture 
has double twist in the nematic field as found in blue phases of cholesteric liquid crystals 
[22,23]. The integrated trajectories are also known to be the projections of great circles of 
S3, the surface of a 4-dimensional sphere, onto flat 3-dimensional space. The curvature of 
S3 implies that the coordination number is less than z = 6 [23,24] which is consistent with 
(4.9). This is illustrated in Figure 3. Here the molecular centers (represented by dots) 
lie on a local triangular lattice, as in the hexatic phase. When the polymers are tipped 
over, the distance between their centers is necessarily increased. This reduces the allowed 
packing density leading to coordination numbers less than z = 6. 

5. The Equivalence of Smectics and Chiral Hexatics 

5.1. Derivation of Chiral Hexatic Free Energy 

I have shown that for K2 > K3 a uniform state with a smectic-C* texture can persist 
for sufficiently strong hexatic ordering. A different scenario is also possible, leading to a 
defect state, analogous to the Abrikosov flux lattice in superconductors and isomorphic to 
the twist grain boundary phase of chiral smectic-A [6]. Consider first the Landau theory 
for an isotropic hexatic in a background nematic field 

^ hex = ^|D^ 6 | 2 + r|^ 6 | 2 + n|Ve| 4 (5.1) 

where ipQ is the hexatic order parameter, which, in the broken phase, is approximately 
ipe ks (v / ^4/6)e 6l6>6 , and, as usual, r oc (T — T c ). D refers to the covariant derivative, 
which subtracts O from the gradient of 9q. 
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Adding to this the chiral term (covariantly modified) in J r g e , I have 

36 

= 6ig fi • [^D^e - VeD^] + y% 2 |^| 2 (A 2 - l) (5.2) 

Recall that under n — > — n, -i/^ <-> V'S an< ^ so the fi rs t term in (5.2) respects the nematic 
symmetry. The second term, though identically zero, will be split up so that I may write 
the total free energy as 

Fhex = J d 3 x |^|(D-6zQ n)^ 6 | 2 + [r-18q 2 ] |^ 6 | 2 + w|V> 6 | 4 + ^nj (5.3) 

Thus, after shifting T c upwards by 18<7g, the free energy of a chiral hexatic is precisely 
that of a chiral smectic [25] where 27r/6(/o, the spacing between surfaces with the same 
bond-order orientation, is the equilibrium smectic layer spacing. There is an important 
difference between the smectic and the chiral hexatic which is what allowed the non- 
chiral ferroelectric term (3.10) in the latter. In a smectic, the physics is invariant under 
Qo —Qo an d 0q — > —6*6 (V>6 V'eO- While that appears to be a symmetry of (5.3) recall 
that 06 — > —9q must be accompanied by n — > — n. In other words, there is a difference 
between qo and —qo in the hexatic that does not exist in the analogous smectic. It is also 
important to note the distinction between O and n. Shifts of 9q by u accompanied by shifts 
of O by du are true gauge transformations. The invariance under these transformations 
is not the sign of a symmetry, but rather a redundancy of identical physical descriptions. 
Rigid rotations in space, on the other hand, are responsible for the gauge-like coupling 
of 6> 6 to n. Under rotation by angle <&, n^n + $xn and 6> 6 — > 6> 6 + [$ x n] • r, so 
that, to lowest order in fluctuations of n around a spatially-uniform, ordered state (e.g., 
n = £), the rotation is a symmetry of the theory. A rotation is physical and corresponds to 
different but equivalent physical systems [26]. In the highly aligned limit, when 5n is small 
and n z w 1, V x fl = 0(5n 2 ) and hence the contribution from the covariant derivative is 
much smaller than that from the chiral term in the derivative %5n. Thus I expect that the 
energetics of the chiral hexatic phases should be similar to that of the analogous smectic 
phases. 

It is now clear that the uniform conical state is just one possible smectic phase. In 
the nematic state (sin0 = 1), there is pure chiral hexatic order, with 6> 6 = q z. This is 
equivalent to a smectic-A phase as shown in Figure 4. Upon reduction of the nematic order, 
the nematic director tips out of the z direction and the phase is equivalent to a smectic-C* 
phase. The chiral order acquires a longer pitch, which is equivalent to a larger smectic 
layer spacing. Finally, as the hexatic order is reduced further, the director lies down into 
a nematic texture and there is no chiral hexatic order. This is just a pure cholesteric. The 
smectic layer spacing has gone to infinity and there is no smectic order. 
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5.2. Defect Phase: The Renn-Lubensky State 

Ignoring the complications of the Mermin-Ho relation and the covariant derivative, 
this model has been studied extensively [6] and, in the type-II regime (where u is sufficiently 
large) one expects a ground state with a proliferation of defects. I propose then a new twist 
grain boundary state of hexatics. It will consist of regions of twisted N+6 separated by 
grain boundaries made up of hexatic disclinations. Within in each region the bond order 
will twist with pitch 2n /% and the director will be well aligned in a nematic state. As one 
moves along a pitch axis perpendicular to the nematic director, across a grain boundary, 
the nematic direction will jump by some finite angle a. This will lead to a state which, 
at long distances, will appear to be pure cholesteric, but will in fact have regions with 
rotating hexatic order. 

Since the possibility of the defect phase is independent of the existence of the uniform 
conical phase, it is possible that the uniform conical phase could be punctuated itself by 
defects. In this case (j)^ would be the equilibrium conical angle and the director would 
relax to this cone. This exotic phase would be similar to the TGBc phase [27] in which 
the defect-free regions are not smectic-A, but rather smectic-C, though in this case it may 
be more appropriate to think of the clean regions as being smectic-C* instead. 

One might speculate on other possible phases. Since the bond order field 6> 6 has much 
the same behavior as $2, a biaxial order parameter, structures of chiral N+6 phases should 
be in close analogy with the myriad of blue phase textures: in other words an N+6 phase 
is closely related to a biaxial nematic. Owing to the connection with smectics it is a 
savory proposition to think of blue phases made of smectic layers, where now the phase 62 
present in the high-chirality limit of a blue phase [22] could be interpreted as ticking off 
the lamellar layers of the smectic in space. This prospect is under investigation [28] . 

5.3. The Structure of Topological Defects 

In light of the Mermin-Ho relation, the structure of the defects requires some clarifi- 
cation. Configurations which satisfy the Mermin-Ho relation can have disclinations in 9q 
without a logarithmically divergent energy per unit length. However, configurations are 
not required to be so innocuous. The defects necessary for a TGB state are still possible. 
The energy of the core will be modified, of course, by the extra energy associated with the 
"curvature" of the nematic field. Considering a screw dislocation in the hexatic order, 9 6 
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will wind around in each xy-plane, independent of z, and so 9q = 6, where 9 is the az- 
imuthal coordinate. The equations of motion imply that at radii larger than a/ K\/Ka/ Qo 
the nematic field will lock to the "superfluid velocity" v = V ± 9 6 /q , and thus 

5n « —9 (5.4) 
q r 

where 9 = (xy — yx) jr is the unit azimuthal vector in cylindrical coordinates. Note that 
since 5n does not depend on z the only component of the Mermin-Ho relation that is 
non-zero is the z component. Using (4.9) I find, to leading order in Sn, that 

[VxOL = -ii (5.5) 

Thus for large r the gauge field O falls off like r~ 3 , which is much faster than the rate at 
which the vorticity falls off (r _1 ) and equal to the rate at which the difference (V ± 9e,— qoSn) 
falls off at infinity. Thus, ignoring the effect of the nematic curvature adds an energy 
comparable to the energy of the defect itself, and certainly does not cause the defect to 
become energetically prohibited. 

Turning back to the Meyer ansatz, one might consider defects in the uniform conical 
state. I can solve the Mermin-Ho relation to find a family of defects in a conical state. 
Using a radially dependent ansatz, I take 

n = [cos (p(r) cos q z, cos <j)(r) sin q z, sin 4>{r)] (5.6) 

In the center of the defect, 0(0) = and I have a pure cholesteric region - a region which 
does not support hexatic order. Moreover, at infinity, the director relaxes to a conical 
state, i.e. 0(oo) = 0<x>- In this case, it is straightforward to compute 

yq cos (p d(p 
[Vxv] x = — - 

= ^gocos0g0 
y r Q r 

r „ , COS 90 

solving for v I have 

sin (y ( t* ) ^ 

v = + +q (sin ^ - sin <f>(r)) z + V9 Q (5.8) 

r 
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where 6 = (xy — yx)/r is the unit azimuthal vector in cylindrical coordinates, and V6> 6 
is a non-singular variation of the bond-angle direction, with equilibrium value gosin^ooi 
in a pure conical state. Far from the core, as r — > oo, v — > V6>6 = qosincfrooZ and is not 
singular. Inside the core, sin^ and, minimizing the free energy for V6>6, the bond 
order is constant. To be more precise, in order to have nonsingular field configurations, 
V x v must not diverge at the origin. Since cos0(r — > 0) = 1, (5.7) implies that d r <p must 
go to no slower than linearly in r. Hence sin <j>(r)/r is well behaved at r = 0. In the core 
of this defect there is pure cholesteric order of the nematic field, while outside the core 
there is a pure uniform conical state. 

In considering configurations (5.8), it is essential to note that, when integrated around 
a closed loop of radius R in an arbitrary constant- z slice, 9q must be single-valued up to 
shifts by 2it/6 and hence 

2im f f 

— = d£-v= dxdy [V x v] 2 = 2tt [sin <j>(R) -sin 0(0)] (5.9) 
" JdM Jm 

Thus, if I consider configurations without any free disclinations, then sin^(r) must take 
on values which are integer multiples of 1/6. This would prevent n from relaxing back to 
a nematic configuration on a longer length scale than the disappearance of In other 
words, the quantization of the circulation of 9q enforces a quantization in the tilt of n. 
Therefore, if there were a core in which ipg vanished (of radius ~ £), outside that core the 
nematic would have to have a constant tilt. This is as if the London penetration depth A 
were smaller than £. This prevents the usual type II defects of the TGB state. This sort 
of restriction is reminiscent of the disclination buckling transition in hexatic membranes 
[29]. This type of defect has a core in which there is not hexatic order, or perhaps an 
annular region in which the tilt of the nematic changes gradually by 27rm/6, for integer m. 
Outside the core there is hexatic order as well as nematic order. Inside there is a gain in 
cholesteric energy as the nematic tips over and exploits its chirality, while there is a loss of 
hexatic energy as its order is destroyed. A detailed energetic calculation, using the entire 
rotationally invariant theory must be performed in order to determine whether this non 
topological defect is stable. 
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Figure Captions 

Fig. 1. Escape into the third dimension of a biaxial object. As a +1 nematic disclination 
tries to escape into the third dimension, the small axis director is forced to have 
a +1 disclination. The long thicker line is the nematic axis while the shorter, 
thinner line is the small director. At the boundary the nematic is normal and the 
small director is uniform and normal to the plane of the boundary. 

Fig. 2. A braided state with a double twist texture. The center line is straight and 
has fewer than z = 6 nearest neighbors. The lines may represent long polymers 
or merely curves with tangent vectors equal to the local nematic director (the 
integral curves). 

Fig. 3. Illustration of the frustrated packing of tilted molecules. Note that when the 
molecules are tipped over, the distance between their centers increases. This 
reduces the allowed packing density and leads to coordination numbers less than 
z = 6. 

Fig. 4. Model of a chiral hexatic. In each plane the bond-order parameter is 6> 6 = 
Oq mod 2tt/6. Between the planes the bond order uniformly precesses along 
the average nematic director, n = z. The planes are analogous to smectic planes, 
though there is no density wave in this liquid crystalline phase. 
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